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Abstract. In this paper, by using two different techniques we derive an explicit formula for the mean 
first-passage time (MFPT) between any pair of nodes on a general undirected network, which is expressed 
in terms of eigenvalues and eigenvectors of an associated matrix similar to the transition matrix. We then 
apply the formula to derive a lower bound for the MFPT to arrive at a given node with the starting point 
chosen from the stationary distribution over the set of nodes. We show that for a correlated scale-free 
network of size N with a degree distribution P(d) ~ d -7 , the scaling of the lower bound is JV 1_1//7 . Also, 
we provide a simple derivation for an eigentime identity. Our work leads to a comprehensive understanding 
of recent results about random walks on complex networks, especially on scale-free networks. 



PACS. 05.40. Fb Random walks and Levy nights 
Classical transport 
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1 Introduction 

As a new interdisciplinary subject field [I], complex net- 
works received considerable interest from the scientific 
community in the past decade [2]|3]|i]. A central issue in 
this area is to uncover the effects of structural properties 
on dynamical processes on networks [MS]- As a paradig- 
matic dynamical process, random walks on complex net- 
works have been investigated intensively {7lt8ll9i rT0lfTTirT2l 
[PSlfTl] . In particular, considerable attention has been paid 
to evaluating the mean first-passage time (MFPT) from 
one node to another in networks, which is required in var- 
ious settings [TS| and has been simulated by studies of re- 
lated problems on complex networks such as transport (TT] 
and search [IMT7JTB]. 

The MFPT TJy from node i to node j is the expected 
time taken by a walker to first reach node j starting from 
node i. Recently, there has been increasing interest in de- 
termining the MFPT (Tj) to a given node j, defined as the 
average of T\j over all starting points i, and in finding the 
dependence of (Tj) on network size N , since this quantity 
is widely used to characterize transport efficiency [T9,20, 
|H p 2 i r ^ |M l |2 M 2 l |2 ^^ . Previous 

works have unveiled nontrivial scaling of MFPT (Tj) in 
a network with size N. For example, it was shown that 
in some scale-free networks the MFPT (Tj) to a most- 
connected node j behaves sublinearly with the network 
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size N [26 ( 27 ( 28 , which is in sharp contrast to that in 
standard regular fractals (e.g., Sierpinski gaskets [2"D1I2"T1 
[22] and T- fractals I23 ( [2"4" [ [2~5j. where (Tj) scales superlin- 
early with N. This striking finding has inspired a growing 
number of works with an attempt to explore the impacts of 
other topological features on the scaling of MFPT (Tj) to 
a target node j, including modularity [33] . fractality |341 
133] . and so forth. However, various topological properties 
are closely related to each other, e.g., scale-free networks 
are often simultaneously small- world [36] , and fractal char- 
acteristic is frequently accompanied by a modular struc- 
ture [37] ■ K is thus very hard to distinguish which indi- 
vidual property plays a crucial role in determining the 
behavior of MFPT (Tj) in a network. In some cases, it 
is possible to confuse and even misinterpret the effects of 
different topological features. 

In this paper, we present a general framework for study- 
ing discrete-time random walks on networks, which can 
unify previous results about MFPT. For this purpose, we 
derive an exact formula of MFPT between two nodes on 
an arbitrary network, in terms of eigenvalues and eigen- 
vectors of a matrix similar to the transition probability 
matrix associated with a random walk. Based on this for- 
mula, we further provide a lower bound of MFPT to a 
target for a walker, which started from a node selected 
from the stationary distribution over all nodes in the net- 
work. We find that the lower bound is proportional to 
the reciprocal of the degree of the target node, regardless 
of other topological properties. Particularly, in correlated 
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networks with N nodes obeying a power-law degree dis- 
tribution P(d) ~ d~ J , the dominating behavior of the low 
bound of MFPT to a hub node is N 1 " 1 ^. To that end, 
we make a comprehensive analysis of existing results about 
MFPT in different networks and show that those results 
can be understood within this framework. Finally, as a 
byproduct, we derive Kemeny's constant |38| for random 
walks on networks. 



where / denotes the identity matrix of order TV, and W is 
the matrix with each row being tt t . Specifically [41] 



(2) 



Next, using the spectral graph theory [32], we show 
that Tij can be expressed in terms of the eigenvalues and 
eigenvectors of a matrix related to the random walks. 
Since P is asymmetric except when graph G is regular, 
we introduce the following matrix: 



2 Brief introduction to random walks on 
networks 

We study discrete-time unbiased random walks [7] on an 
arbitrary connected undirected network G with N nodes 
and M edges. The N nodes of G are denoted by 1, 2, 3, . . . , N 
The connectivity of G is encoded in its adjacency matrix 
A, whose entry a,ij = 1 (or 0) if nodes i and j are (not) con- 



N 



S = D-iAD-i = D*PD~2. (3) 

Obviously, S is real and symmetric and has the same set of 
eigenvalues as P. Moreover, if ip is an eigenvector of S cor- 
responding to eigenvalue A, then D~2ip is an eigenvector 
of P associated with eigenvalue A. 

Let Ai, A2, A3, • • • , Xjsr be the N eigenvalues of matrix 
S, rearranged as 1 = Ai > A 2 > A 3 > . . . > \jy > — 1, 
and let ipi, ip2, ip3, ■ ■ ■, ipN be the corresponding mutu- 
ally orthogonal eigenvectors of unit length, where ipi = 
(ipn,ipi2, ■ ■ .,ipiN) T - Then 

\P T S\P = diag[Ai, A 2 , 



• A 



\p\p T = \p T \p = I 



Thus, the entry ipy of & satisfy 



TV 



nected by an edge. The degree of node i is di = 
and the diagonal degree matrix of G, denoted by D, is 
defined as: the ith diagonal element is di, while all non- 
diagonal elements are zero. Thus, the total degree of all 

nodes is K = 2M = J2f=i di, and the average degree is where & is an orthogonal matrix, that is, 
(d) = 2M/N. 

In the process of a random walk considered here, at 
each time step the walker starting from its current location 
moves with a uniform probability to one of its neighboring 
nodes. In fact, this stochastic process of random walks is 
characterized by the transition matrix P = D~ 1 A, whose 
element pij = a^jdi presents the probability of jumping 
from i to j in one step. Let P express the tth power 
of matrix P. Its ijih entry denoted by (p l )ij represents 
the probability for a walker, starting from node i, to visit 
node j in t steps. It is easy to verify that such a random 
walk is actually an ergodic Markov chain [31?ll4"0] . whose 
stationary distribution tt = (tti,tt2, ■ ■ ■ , ttn) T is a unique 
probabilistic vector satisfying m = di/K, J2 i=1 ^i = 1, 
and tt t P = 7r T , where the superscript T means transpose. 



N 



1) if i = 3, 
0, otherwise. 



(4) 
(5) 

(0) 



k=l k=l 

Equation (U) can be rewritten in an alternative way as 

5 = ^diag[A 1 ,A 2 ,...,A A ,]<f T , (7) 



which indicates that the entry sy of matrix S has the 
following spectral representation: 



N 



Sij = AfcVfeiV'. 



kj ■ 



(8) 



k=l 



3 MFPT between two nodes 

A quantity of primary interest related to random walks is 
MFPT. It encodes a great deal of information about the 
random- walk dynamics [TS]- In what follows, we use two 
different approaches to derive MFPT between an arbitrary 
pair of nodes. 



3.1 Method of spectral graph theory 



It is well known that the MFPT from node i to node 
j can be expressed exactly in terms of the elements of 
the fundamental matrix Z for the corresponding Markov 
process, defined by 



We proceed to express matrices /, P, and W in terms 
of ^ and D. First, according to Eq. ([5]), one has 

/ = Z)-2^diag[l,l,...,l]tf rT Z)2. (9) 

On the other hand, from Eqs. ([3]) and ([5]) one obtains 

P = D~iSDi = D^iFdiag [A x , A 2 , . . . , X N ] <P T . 

Finally, it follows from the definition of the stationary 
probability [591I4H] that 

W = lim P t = lim ( D~$SD^ * 



z= (i -p + wy\ 



a) 



lim£>-i [Va^jW 
\fc=l / 

D-i^jDi +)imD-i ^A^fe^j V h 

D-3^diag[l,0,...,0]^ T Z>3 5 (n) 
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where we have used the facts Ai = 1 and |A*| < 1 (k > 2) 
when G is not a bipartite graph. 

The above three equations lead to 



Z = D 2 i^diag 



1. 



1 



1 



1-A 2 

Then, rewrite the entry Zij as 



1 - A 



N 



<P f Dz. (12) 



ij 1 



% k=2 \ K 



w dU 



(13) 



In order to express MFPT in terms of the eigen- 
values and eigenvectors of matrix S, all that is left is to 
determine the eigenvector tpi of S corresponding to eigen- 
value Ai = 1. Since S = D^PD^? and Sipi = ip%, one 
has P(D~2ipi) = {D~2ipi). Since PI = 1, one also has 
_D~2-0i = 1, so that 



K 



(14) 



Inserting Eqs. (TT3l and (Ti"4"f into Eq. @ and with some 
calculation one obtains 



Ti 



-*f-L-(£ 



(15) 



which will be very useful in the following derivation of 
MFPT (Tj) to a given target node j. 



3.2 Method of generating functions 



Notice that Eq. (|T5|) can also be derived by using the for- 
malism of generating functions [33]. Let {<f)ij represent 
the probability of a walker from node i to first reach node 
j after t steps. Note that (<?*)ij is quite different from the 
(j? )ij studied above, since the latter represents the prob- 
ability of hitting j in t steps given that the walker stated 
from node i. The two quantities satisfy the following rela- 
tion: 



(16) 



s=0 



This relation can also be expressed in terms of generating 
functions. To show this, we define 



and 



Qij{x) = ^{q^ijZ* , 
t=o 



(17) 



(18) 



Pij{x) < 1, while Qij(x) is continuous for x G [0,1] but 
may diverge as x — > l~. Equations (fTF]) and (fT8]) together 
yield a useful relation given by 



(19) 



By its definition, the MFPT Tij can be obtained via 
differentiating Qij{x) with respect to x and then setting 
x = l: 

d 



T — 



(20) 



Thus, our goal is reduced to proving that Eq. (|20j) is equiv- 
alent to Eq. (fl5|) . 

As shown in Eq. (fTTj) . one has 



so 



AT 



0%- = £ A fcV , fe^. 



fc=l 



(21) 



(22) 



Thus, for < 1, one has |Afca;| < 1 (k > 2) and 

oo Af 



/=o fc=i 

TV 
fe=l 



di 

N oo 



EE( Afca; )*^ w ^ 



k=l t=l 
Af 



^WxE^+E 



t=l fc=2 
Af 



'•7 1 



\ - AfcX /gL- 

-j h V -j \ — WkiVkj\ -r, 

1 — x 1 — Xhx \ di 



(23) 



where we have used Eq. (0 and i/'iiV' 
that 

XkX 1 1 — AfcX 



WV * = 7ri ' Note 



1 - A fc a; 1 - A fe .x 1-Afcar' 

the second term, denoted by ^2, on the last line of Eq. 
can be recast as 



(24) 



.92 



N l /cF w l 

fc=2 K " 1 fc=2 



fe=2 
Af 



i - A fc x /d7 

t — tpuwkj\ -r 

X k x V di 



fc=2 
Af 



k=2 

Therefore, 



(25) 



where a: is a complex number. By the Abel theorem, Py (x) 
is a uniformly continuous function for a; G [0, 1] and < 



1 * 1 

fc=2 



^7 



(26) 
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Analogously, Pjj(x) can be derived. Applying Eq. (|22l) Then the derivative of Qij(x) with respect to variable x 



gives 



oo N 



t=0 k=l 
N 



N 



*!=i 



fc=l t=l 



N 



fc=i 



N 



1+7T,-- + 

1 — X 



St 

fc=2 



fc=2 



AfcX o 
"V'fc 



(27) 



is At x = 1, all items with (a; — 1) will be zero. So, 



JV 



■"■i E rr^^fe v | + ^ E i^-'t 



fc=2 



fc=2 



3 k=2 K \ 

1 - A fc I eL 

fc=2 K \ J 



(32) 



which is consistent with Eq. (| 15[) . Notice that the formula 
The last term, denoted by g 3 , in the last line of Eq. (23 in EqS " ^ and ® haS ^ previousl y derived 



can be further written as 



in [44] by using a different approach. 



N 1 * 1 

ft» = Ei3jri^-Ei 



fc=2 
JV 



fc=2 



1 - X k X 2 

~\k~x i 



Y — 



k=1 
N 



Y^— 



1 + 

1 + TTj 



fc=2 



Plugging Eq. flSS) into Eq. (J27J) leads to 



1 ^ 1 



fe=2 



4 MFPT to a given node 



After deriving the MFPT Ty from node i to node j in a 
general network G, we now study the trapping problem 
defined on G, a special random walk with a deep trap lo- 
cated at an arbitrarily given node, e.g., node j. Let (Tj) 

(28) 

be the average trapping time (ATT), which is the mean 
of MFPT Tij to the trap node j, taken over all starting 
points. There are two slightly different definitions for (Tj): 
For the first one, the average (•) is taken over the uniform 
distribution of the set of starting points [19112011211 1521125 ] ■ 
while for the other, the average (•) is taken over the sta- 
tionary distribution |12] . Since both definitions lead to an 
identical scaling of (Tj ) with network size N [T2J , here we 
(29) adopt the second one, given by 



Substituting the above obtained expressions of Pij(x) 
and Pjj (x) into Eq. p^|) yields 



Qtj yxj 



N 1 

k=2 _ 

N 1 

TTj-r— + y 

1 - , 



k=2 
N 



- Xk.x 



^■+( i -^Et^^^V| 

fc=2 ' 

N 1 



(30) 



f N 

(Tj) = ~ _ ~ E ^ T n ■ 



(33) 



Substituting the expression of in Eq. (TTSJ) to Eq. (|3"51) 

gives 



(Tj) = 



i N N i (4? 

i n a n i / 

1 - TTj f-f dj 1 - Afe I 

i=l fc=2 K V 

1 N ( 1 W J 

1 \ "« / t 2 \ ^ "i 

-—£(— 

1 - 7T, ^ I 1 - , 

3 k=2 \ 



4>kj - i>kri>k 3 \ -j- 



1 V *' 

-V'fcjWxE^V^ ' ^ 

V J i=l * 
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N i N I "I r N 



i 

^(^Et^t^ 



" N 1 * A 

fc=2 fc=2 v K ' 




JV 1 r^~" 

2 V Gta 




~ 1 

fc=2 


2 



(31) 



It then follows from Eq. © that Yh=i Vfe = J2iLi V'fe-. 
Thus, the second term is equal to zero. So 

= 1 1 ,2 

l-TT^.^l-Afc^' 



(35) 



We next derive a lower bound for (Tj). By Cauthy's 
inequality, one has 

By Eq. ([5]). in matrix S 1 , each of its diagonal entry Sjj = 

T,k=i X k^lj > (J = 1,2,..., TV). On the other hand, 
the traces of matrices S and P are both equal to zero, so 
Sjj = 0. Thus, since Ai = 1, one has 



N 



N 



N 



2(i - A fe )vi- - EC 1 - = 1 - E = 1 • 

(37) 



fc=2 

In addition 



fc=i 



fc=i 



N 



JV 



E ^ = E - = i - ^ 



fe=2 



fc=l 



Therefore, 



and 



fe=2 



1 K.. , 2 K. 



K 



(38) 



(39) 



(40) 



This lower bound is sharp in the sense that it can be 
achieved in some graphs, e.g., a complete graph with JV 
nodes where the MFPT from an arbitrary node to node j 
is JV — 1 [pj , which is exactly equal to jf- — 1 since in this 

case K = N(N — 1) and dj = N — 1. 



^ li= QNote that most real networks are sparse with the aver- 
age degree (d) — ^ being a small constant. Thus, in these 
networks, for a target node with degree d, the lower bound 
becomes N ^ — 1, which, in the limit of large network size 
JV, is proportional to JV and the inverse degree of the trap 
node. Moreover, it has been realized that most real-world 
networks display ubiquitous degree correlations [45ll46j . At 
the same time, they are scale-free with their degree distri- 
butions P(d) following a power-law form P(d) ~ <i -7 
3,4 , implying that there exists hub nodes in these cor- 
related networks whose degree is 

[6] . These nodes play a dominant role in various dynamics 
running on networks [5,6 , including the trapping prob- 
lem addressed here. Noticeably, when a trap is fixed at 
a hub node with degree d max , Eq. (|4"0|) shows that the 
lower bound for ATT is — 1. In large networks (i.e., 

JV -> oo) the scaling of ATT is JV 1 " 1 /^" 1 ), which grows 
sublinearly with the network size JV. 

To this end, we have presented a unifying framework to 
encompass many existing works about ATT on networks, 
especially on scale-free networks with degree correlations 
among nodes. First, in a generic scale-free network, al- 
though the low bound [Eq. (|3U|)] for ATT to a hub node 
is hard to reach, but the minimal scaling JV 1-1 /^ 7-1 ) can 
be achieved in some particular scale-free networks, e.g., 
hierarchical and modular scale-free networks [321133] that 
thus have the most efficient structure for transport by dif- 
fusion when the target is positioned on a hub, even on a 
particular non-hub node [15] ■ Second, in some scale-free 
networks, although the scaling for ATT to a hub increases 
sublinearly with the network size, it cannot achieve the 
minimal scaling JV 1- VOf - !■). Previously studied network 
examples include the pseudofractal web [37] , the Apollo- 
nian network [25] , and the (1,3)— flower [30]. Third, the 
ATT to a hub can also display a linear scaling in some 
scale- free networks, such as the Koch networks [29 1 131 ] and 
the (2,2)— flower [30]. Finally, in a certain scale-free net- 
work family, the ATT to a hub can even scale superlinearly 
with the network size [3~4"ll35| . 



5 Alternative derivation of eigentime identity 



As can be easily seen from Eq. ([55]). the MFPT (Tj) to 
a target node j is strongly affected by the degree of the 
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target node j. However, the average of MFPT (Ty) from 
node i to a node j randomly chosen from all nodes ac- 
cordingly to the stationary distribution is independent of 
the starting point i. Let (T,-) denote this quantity corre- 
sponding to starting node i. Then 

N 

(T i ) = '£w j T ij . (41) 
i=i 

It was shown in [3D] that the sum 53i=i n j ^ij does not 
depend on i. Using different methods [38 , 40 ,49 , it has 
been shown that (T<) is a constant given by 

^> = Ei^' (42) 

which is referred to as an eigentime identity [10J, some- 
times Kemeny's constant |38j . 

Below we use an alternative yet simple way to re-derive 
Eq. (J42J) . According to Eqs. © and ([T5]). one has 

N N 

( T i) = ^2^j T v = E(%' - 




It then follows from Eq. © that J2j=i ^/y = 1 anc ^ Hj=i ^< 
E^Li i>kji/)ij=0- Thus, Eq. (O reduces to Eq. (gSJ). 

6 Conclusions 

In this paper, we studied theoretically some aspects of ran- 
dom walks on general networks. A discrete-time random 
walk on a network is described by its transmission matrix, 
and thus information about the random walk process is 
encoded in some related matrices. Through applying two 
different approaches, we demonstrated how the MFPT be- 
tween two nodes can be expressed in terms of the eigen- 
values and eigenvectors of a matrix similar to the trans- 
mission matrix. Then, we derived an exact formula for 
MFPT to a given target node by expressing it in terms of 
the eigenvalues and eigenvectors of the same matrix, based 
on which we further provided a low bound for this inter- 
esting quantity. We showed that the low bound depends 
only on the degree of the target node and is seemingly 
independent of the structural properties of the network. 
Particularly, for scale-free networks with degree correla- 
tions, we expressed the minimal scaling of this low bound 



for a most-connected target node by the degree distribu- 
tion exponent 7, which displays a sublinear dependence 
on the network size. This provides a broader view of pre- 
viously reported results for MFPT to a given node on 
different networks, especially on scale-free networks. 

In addition to the above-mentioned main results, some 
of our byproduct results are also very useful. For exam- 
ple, we gave a spectral representation for the entries of the 
fundamental matrix of the studied random walks. Making 
use of this expression, we provided a simple derivation 
for the eigentime identity. As another application, this 
spectral representation can be used to prove the follow- 
ing symmetry property about random walks on a network: 
Tab + T bc + T ca = T ac + T cb + T ba , which holds for three ar- 
bitrary nodes (a, b, and c) and was first discovered in |50) . 
We expect that this important result can find many other 
applications. Finally, we believe that this work can help 
better understand random walks on general networks, to 
provide useful insights in related dynamical processes on 
the networks. 
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